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ABSTRACT 

Three-parameter eigenvalue problems are discussed in this paper. Tensor Rayleigh Quotient Iterative Method for 
Three-Parameter eigenvalue problems are discussed in this paper. Finally some numerical results are presented to illustrate 
the performance and application of this method. 
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1.1 INTRODUCTION 

Multiparameter eigenvalue problems are generalization of one-parameter eigenvalue problems and can be found 
when the method of separation of variables is applied to certain boundary value problems associated with partial 
differential equations. Much more works have been done in the field of one -parameter eigenvalue problems, both 
theoretically and numerically compared to two-parameter or more than two-parameter eigenvalue problems. Some works 
have been done theoretically in the field of multiparameter eigenvalue problems. Few authors have dealt with the 
multiparameter eigenvalue problems numerically mainly in two-parametic cases. Numerical methods applied to a three- 
parameter problems are very limited and hence some contribution in this area are always in needed. 

1.2 THREE-PARAMETER EIGENVALUE PROBLEM AND ITS REDUCTION TO A SYSTEM OF 
ONE-PARAMETER PROBLEMS 

Consider the three-parameter eigenvalue problems 


A w x = AjA n x + A 2 A l2 x + Ay A u x 

^20 y — y ^2^22 y 

A 3Q z = A i A 31 z + A 2 A 32 z + A 3 A 33 z (1.2.1) 


Where A- 6 , i= 1 ,2,3 and 
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xe [ 1 ;- ! \{0}, A 10 ,A n ,A 12 , A 13 e □ "*" 

0},A 20 ,A 2I ,A 22 ,A 23 ea" 
ze □ ^ \{0}, A 30 , A 31 , A 32 , A 33 € □ pxp 

Where A, E H , i= 1 ,2,3 are called the eigenvalues and x,y,z are called eigenvectors of the problem. 
Problem (1.2.1) can be reduced to a system of three one -parameter problems: 


A X U = /IjAqM 
A 2 u = A 2 A 0 u 
A 3 w = A 3 A 0 u 

where A 0 , A p A,, A 3 are (mnp)X (mnp) dimensional matrices defined ; 


( 1 . 2 . 2 ) 


A) At®A2®As A 1 0 As 0 A 2 + A 2 0 As 0 At A2 ® At ® A: 


-*"^13®Al ® A 2 As 0 A 2 0 Al 


A ^ 10 ® A2 ® A 3 Ao ® A 3 ® A2 A2 ® A3 ® Ao A2 ® Ao ® A: 


+A ]3 ® A, 0 ® A 32 A 13 ® A^ ® A 30 


a 2 = a 11 ®a w ® a 33 - a u ® A 3 ® Ao + Ao ® As ® At - Ao ® At ® A: 


+A 3 ® At ® Ao A 3 ® Ao ® A 


A 3 A u ®A 22 ®A 30 Ai 0 Ao 0 A2 A2 0 Ao 0 At A 12 ®A 21 ®A 3( 
+Ao ® At ® A2 — Ao ® A2 ® At 


u - x09 y 09 z 


Theorem : Let (A, A 5 A) be an eigenvalue and (x, y, z) a corresponding eigenvector of the system (1.2.1) then 
(A, A, A) is an eigenvalue of the system (1.2.2) and u — X® y ® Z is the corresponding eigenvector. 


Definition 1.3.1. The Kronecker product (. ® .) : □ mx " XD p q — > Q mpx " q j s defined by 
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a n B ■■■ a ln B 


A®B 


a.B ••• aB 

n 1 nn 


Where we use the standard notation (A)y = 

The Kronecker product is a special case of the tensor product, and as such it inherits the properties of bilinearity 
and associativity, i.e. 

(kA) ® B = A® (kB) = k(A® B) 

A® (B+ C) = A® B+ A® C 
(A+ B) ® C = A® C + B®C 

We now establish a famous property of the Kronecker product, from [9]. 

Lemma (Mixed product property). Let Ae □ mxn ,Be □ pxq ,Ce □ nXk ,De □ qXr . Then 
(A ® B){ C® D) = (AC® BD) . 

In particular, if A, B E U "" and X, y E U then 
(A ® B)(x® y) = Ax® By . 

2.1 TENSOR RAYLEIGH QUOTIENT 

The Tensor Rayleigh Quotient p(x, y, Z, A J0 , A n , A 12 , A ]3 , Aj,,, A, 1; A, 9 , A^, A 30 , A 3| , A 3 -,, A 33 ) is an triplet 
(P \ , p 2 , P 3 ) such that 

u t A,u 

A=^r- 

u A 0 u 

u t Am 

Pt- ~ T a 

u A 0 u 

u T Am 
u A 0 u 

Where U = X® y ® Z 
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Hence the Tensor Rayleigh Quotient at an exact eigenvector is 


_ U \u _ \u T \u _ 3 

A — t A — T7 — A 

u A 0 u u A 0 u 


_ u 7 A 2 u _ A^u 7 A 0 u _ 3 

P2 T 

W A 0 w w A 0 m 


_ u r A 3 u _ A^u 7 A 0 u _ 3 

A — r a — 77 — ^3 

u A 0 u u A 0 u 


Now 


u r A 0 u = u t (A u ® A 22 ® A 33 - A n ® A 2 3 ® A 32 + A p ® A 23 ® A 31 - A 12 ® A 21 ® A 33 
+ Ap ® A 21 ® A 32 — Ap ® A 22 ® A 31 ) u 


u r A q u = ( x r ®y T ®z )(A, , ® A 22 ® A 3 , - A, , ® A 23 ® A 32 + A 12 ®A 23 ® A 31 - A 12 ® A 21 ® A 33 
+ A 13 ® A 21 ® A 32 - A 13 ® A 22 ® A 31 )(x® y® z) 

= x 1 A n x ® y T A 22 y ® z T A 33 z - x‘ A n x ® y‘ A 23 y ® z‘ A 32 z + x‘ A l2 x ® y‘ A 23 y ® z‘ A 3l z - x T A n x ® y T A 2l y ® z T A 33 z 
+x T A l3 x ® y T A 2i y ® z A 32 z — x‘ A 13 x ® y T A 22 y ® A 31 z) 

But X 7 A n JC , y r A n y etc are scalars. Hence ® becomes normal multiplication. So 

u 7 A 0 u 

= ( x T A n x)(y T A 22 y)(z T A 33 z) ~{x‘ A n x){y J A 23 y)(z' A 32 z) + (x T A n x)(y r A 23 y)(z T A 3l z) - ( x T A l2 x)(y T A 2l y)(z‘ A 33 z) 
+(x r A l3 x)(y T A 2l y)(z‘ A 32 z) — (x T A l3 x)(y T A 22 y)(z T A 3l z) 


Similarly 


u 7 A jM = (x T A 10 x)(/A 22 ;y)( A A 33 z) - (x‘ \ 0 x)(y T A 23 y)(z T A 32 z) + (x 1 A l2 x)(y T A 23 y)(z T A 3Q z) - (x 7 A n x){y‘ A 20 y)(z r A 33 z) 
+{x A i3 x)(y T A 20 y)(z r A 32 z) - (x 1 A l3 x)(y T A 22 y)(z T A 30 z) 

u T A 2 u = (x T A n x)(y T A 2O y)(z T A 33 z)-(x T A ll x)(y T A 23 y)(z T A 30 z)+(x T A l0 x)(y T A 23 y)(z T A 3l z)-(x T A l0 x)(y T A 21 y)(z T A 33 z) 
+(x T A l3 x)(y T A 2l y)(z r A 30 z)-(x T A l3 x)(y T A 20 y)(z T A 3l z) 

u T A 3 u = (x r A ll x)(y T A 22 y)(z T A 30 z)-(x T A n x)(y T A 20 y)(z T A 32 z)+(x T A l2 x)(y T A 20 y)(z T A 3l z)-(x r A l 2 x)(y T A 2l y)(z T A 30 z) 
+{x T \ 0 x){y T A 2l y){z T A 32 z)-(,x r \ 0 x){y T A^yXz 7 A il z) 

2.1 TENSOR RAYLEIGH QUOTIENT ITERATIVE METHOD 

Consider the three parameter eigenvalue problem (1.2.1). In mtrix form it can be written as 
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A l0 x - A n A t x - A I 2 / l 2 x - A l 3 A 3 x 
A 2Q y — A 2 l A t y — A 22 A 2 y — A 23 A 3 y 
Ao^ — A 3 I A,Z — A 32 A 2 Z — A 33 A 3 Z 





Newton’s method applies to the equation F(u)=0 gives 

F\u k ){u k+l -u k ) = -F{u k ) 


Where 

F'{u) 


Ao AAi AA2 AA3 

0 

0 

-Aix 

-A I2 x 

-A n x 

0 

A 

0 

l 

> 

A 

1 

A 

A 

1 

> 

A 

0 

-Aiy 

— Aj 2 y 

— A^y 

0 

0 

Ao — AAi ~ A A2 ~ AA 3 

-Aiz 

~A 32 z 

~A 33 z 

-x T 

2 

0 

0 

0 

0 

0 

0 

A 

0 

0 

0 

0 

0 

0 

1 r 

— z 

0 

0 

0 


2 


Now to apply Tensor Rayleigh Quotient Iterative method for three-parameter eigenvalue problem we replace 
respectively A , A , A with p { ,p 2 ,p 3 in equation (1.2.1). Let (x 0 ,y 0 ,Z Q ) be the eigenvectors where 
||x 0 || = 1, ||y 0 || = 1, ||z 0 1| = 1 for the problem (1.2.1). An initial approximation for the eigenvalue is not needed. Since we 
can calculate one using the Tensor Rayleigh Quotient. 

Now we have to solve 


F(u)=0 where 


Ao x — A iPi x ~ \iPi x — P\\P\ X 
Ao y — AnPj — \ 2 P 2 y- A 22>Ppy 
A 30 z — A 3 x p x z — A n p 2 z — A 33 p 3 z 


F(u)= 


1 r 1 

— x x-l 

2 

1 ' 1 

-V y-1 


1 


1 1 

z k z k -\ 


Using Newton’s method we have 
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F'(u k )(u k+1 -u k ) = - F(u k ) (2.1.1) 

Where 


> 

o 

1 

35 

> 

1 

35 

to 

> 

1 

,35 

> 

0 

0 

~Ai x 

-A l2 x 

-A l3 .r 

0 

A 2 o — P\ A>i — PiAi — PiAii 

0 

-Aiy 

~Ai y 

~ Aiy 

0 

0 

© 

1 

35 

Jfc- 

1 

35 

A 

i 

35 

Jfc- 

— A 3[ z 

~ Aiz 

— A 33 z 

-x T 

2 

0 

0 

0 

0 

0 

0 

\ f 

0 

0 

0 

0 

0 

0 

1 r 

— z 
2 

0 

0 

0 


Thus (2.1.1) => 


Ao -Pi k) Ai -p 2 k) \i - A w A 3 0 0 

o Ao — A Ai ~ Pi A2 — A A3 0 — Ai yk — A 2 ^ — A 3^ 

0 0 Ao “ A Ai - Pi A2 — A A3 — Ai Zk ~A 32 z k —A^Zk 


1 . 

— X 
2 4 

0 

0 


:v t 


0 


1 . 
2 


Ay* 

Az* 

AAf 

AX * 


, (*), 


AoAk — A 1 A At — A 2 A Ak — A 3 A 

Ao } 7 k — AlA yk ~ AlPl yk — A.2P3 yk 

Ao Zk — A I P\ Zk — AlPl Zk — A3P3 ^jfc 


~ X k X k ~ 1 


T 

■y k y k - 



( 2 . 1 . 2 ) 


Consider the 1 st component of (2.1.2) 

(Ao - a w Ai - A K) A 2 - p3 K) A3)(*k +Ax k ) = A A A A (k) + A 2 x APi ] + A 3 ** A A (t) (2>L3) 

Consider the 2 nd component of (2.1.2) 
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(Ao A At Pi A 2 A AsKY*"*” A Ya) AiYa A A A 2 y AP 2 "'‘AsYa^A 


,(K) , 


,(/n 


00 






Consider the 3 rd component of (2.1.2) 


(Ao - /^Al - /?2 K> ^32 - Pl K>A 33 X Z A + Az a ) = Al^VA + A 2 Z A A /A' + A 3 Z A A /A' 


Consider the 4 th component of (2.1.2) 

'*'* X k+\ = 2x k x k — 1 
Similarly from the 5 th and 6 th component 

y T k y k+ x = 2 y k T y k - 1 

T '1 T 1 

Z’k Zk+l Z’k ^ 

From (2.1.3) 

^A + i = (Ao - aA A 1 - aA A 2 - aA A 3 A (A pApV + A 2 *a a aA + A 3 *A A AA) 

Similarly 

Ya+i = ( A 2 0 -pl K)A 2 l -P? ) A 22 -p ( 3 K) A 23 r l (A 2 l y k Apl k) + A 22 y k Ap { 2 k) + A 23 y k ApA) 
Za +1 = (Ao - aAAi - A® A2 - p ( 3 K) A 33 y\A n z k Apl k> + A 32 z k Ap ( 2 k) + A 33 z k Ap ( 3 k) ) 


Multiplying (2.1.9), (2.1. 10), (2.1. ll)on the leftby x\, y \ , z\ and using (2.1.6), (2.1.7), (2.1.8) 


x k T u k 

T 

x k v k 

* 

A 

1 



2x/x,-l 

T 

y k p k 

T 

y k q k 

T 

y k h 

a aA 

= 

2 y k T y k ~ l 

z k T a k 

Zk T b k 

T 

Z k c k 

a aA 


1 

K> 

??• 

1 

1 


Where 

«a = (Ao - aAAi - pfA xl -pf ’As) -1 A a 
v* = (■ Ao - A W Ai - P?A n - A (Jr) A 3 )"‘ A 2 ** 
^ = (Ao - a ot Ai - A K) A 2 - aAAsA \,x k 
Pk = ( Ao - Pi K) At - A K) A 2 - A* 1 A 3 A At ?A 


35 

(2.1.4) 

(2.1.5) 

( 2 . 1 . 6 ) 

(2.1.7) 

( 2 . 1 . 8 ) 

(2.1.9) 

( 2 . 1 . 10 ) 
( 2 . 1 . 11 ) 
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A = ( A 20 - P [ K ' A 21 - P? )A 22 ~ Pr^y' ^k 

r k = (A 20 - p[ K) A 21 - /?fA 22 - p^’A^y' A 23 y k 

A = ( A 30 - A (JC)A 31 - / ? 2 K>A 32 - P^'^sT 1 A nZk 
K = (■ 4o - A^l - />2 K)A 32 - >A 33 )"’ A 3 2 Z* 

A = (4o - Pl K)A 31 ~ py >A 32 - P ( 3 K)A 33T' A 33Z k 

The following algorithm can be used for Tensor Rayleigh Quotient Iterative Method 

• Start with an initial approximations (x Q , y 0 , Z 0 ) to the eigenvectors where ||x 0 1| = 1, || y 0 1| = 1, ||z 0 || = 1 Then 

Calculate the Tensor Rayleigh Quotient (p k , pj , p k ) 

• Check determinant of matrices 

[Ajo — A Al — Pi A 2 — Pi A l 3 ] ’ [ A 20 _ A A 21 ~ Pl A 22 ~ P 3 A 23^i\- A 3Q ~ P \ A 31 — A A 32 — A A 33 ] 
and if either are equal to 0, perturb (p k , p! , p!' ) slightly. 

• Solve the following equations 

[A 10 - A (t) At ~ PT \i~ PV A ^k = A n x, 

[A 10 - A (t) At ~P?\ 2 ~Pi k>A nl v, = A 12 x a , 

[A 10 - A (t) At -pfA 2 - AfAa] w t = A 13 x a 

[A20 -p[ k) A lx -p\ k) A 22 -p (k) A 2 3 ]p A = A 21 y jt 
[A 20 -p[ k) A 2l ^p (k) A 22 -p?% 3 ] q k = A 2 2 y k 

[A20 -p^A.-p^^ -p3 k>A 23K = A 3 .Va- 

[A30 - A W At -pV a 32 - A W A 3 ] A = A 31 z i . 

[A30 — p\ k) A 3 \ -p ik)A 32 - A W A 3 ]A = A 32 Z a 

[A30 - A W At - - A W A 3 ]c, = A 33 z a 


• Set up and solve the following system 


3 

fi 

* 

1 

T 

A v k 

x k w k 

~ Ap[ k) 1 


2.x k T x k -\ 

T 

y k p k 

y k a 

T 

y k r k 

A pV 

= 

2 yJ y k - 1 

e 

hi 

IN? 

1 

T 

y K 

1 

X 


2 z k T z k -i 


• Update the approximate eigenvectors 
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+&Pi K \ +A A k) w k 
y k+ i= A Pl k) p k + A / ? fV i - +^P3 k) n 
z k+l =Api k) a k +A pf\ +A pfc k 

• Normalise the approximated eigenvectors 


x k+ l = 


y k+ 1 = 


z k +i 


X k+ i 

l** + ll 

1 m 

||y * +1 

H+i 
|^jfc+i I 


3.1 NUMERICAK EXAMPLE 


We now present a numerical example to show the behaviour and application of our method 
Consider the three parameter eigenvalue problem 


( 1 0 Vx A 


0 2 


\ X 2j 


= 4 


f 3 OY x ^ 


0 4 


J\ X 2J 


+ Z, 


Y 0 Yx A 


V 0 


+ / 1 , 


oYx, A 


0 8 


J\ X 2 j 


f 2 oY)Y 


10 b 


{y 2 J 




4 8 OYy^ 


v 0 9 j 


\J 2 ) 


+ ^2 


bo 0Y3O 


0 15 




+ /I3 


04 oY )Y 


v 0 13 , 




'3 0 ^ 

to 2, 


f , \ 




■K 


30 0 
0 6 


r „ \ 


V 4 2 y 


+ 


75 0 


v 0 16, 


f „ A 


V 4 2 7 


+ A3 


57 0 

0 14 


f „ \ 


\ C 2 J 


Table 1 


Starting Eigenvector 


Iteration 


^=( P r,Pb\prf 


p iK+i) ~b K) 


x 0 = 


3^0 


-0 = 


f 1 > 

-• 05 , 

r 1 i 

V-- 0 Y 

f 1 1 

v -.05 y 


4 —. 0022 ^ 
-.1501 
2514 


v 


4 -8.9850* -009^ 
-.1500 
.2500 


.0711 
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4 


( 2.07970e-016"* 
0 

v - 142 9 J 


1.4344e-016 

*0 = 
y 0 = 
^0 = 

".02" 

".05" 

v2, 

f 1.2 
,-.0J 

l 


0 


"-.I860" 1 
-.3401 
v .5981 y 




1 


"-.1856^ 
-.3402 
v -5979 , 


4.5826e-004 


2 


"-.1856^ 
-.3402 
v -5979 , 


0 

II II II 

K° ?? 

v2, 

" 1.7 

v--0- 

".02" 

vl-5. 

b 


0 


1.1778" 
-.3929 
v 1.0967 




1 


"-1.5669 
-.5222 
v 1-4115 


.5169 


2 


"-1.7946 
-.5982 
v 1.5956 


.3025 


3 


"-1.8" 
-.6 
V 1-6 y 


.0072 


4 


"-1.8" 
-.6 
v 1-6 , 


0 
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Conclusion: Table 1.3 shows that the successive difference between the eigenvalues are gradually decreases. So 
the method converges to exact solution rapidly. In this method the approximate eigenvalues are obtained easily. So one 
can use this method easily to solve three -parameter eigenvalue problems. Here the approximated eigenvalues are 

(-1. 1924e-016, 1500,. 2500) r , (-.25,-1. 8499 e-016„25) r ,(2.07970e-016,0,.1429) 1 ', (-.1856, -.3402, .5979/, (-1.8, -.6, 1.6/, (-1.6364, -.2727, 1.2727) r , (-.236 1, -.3161,. 6052f 
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